Scarring in open quantum systems 
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We study scarring phenomena in open quantum systems. We show numerical evidence that 
individual resonance eigenstates of an open quantum system present localization around unstable 
short periodic orbits in a similar way as their closed counterparts. The structure of eigenfunctions 
around these classical objects is not destroyed by the opening. This is exposed in a paradigmatic 
system of quantum chaos, the cat map. 
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Open quantum systems are very important in many 
areas of physics. For example, they play a central role 
in the study of quantum to classical correspondence [1], 
microlasers [2, 3], quantum dots [4], chaotic scattering 
[5], and more. However, there are several properties of 
these systems that are less known if compared to those 
of closed ones. 

Quantum evolution in open systems is given by 
nonunitary matrices, whose eigenstates (resonances) are 
nonorthogonal and the eigenvalues are complex with 
modulus less than or equal to one. One of the main con- 
jectures about the properties of the spectrum is that the 
mean density of resonances follows the fractal Weyl law 
[6] . This law predicts that the number of eigenstates that 
have a finite decay rate goes as N 1 cx %~^ x \ where d is 
a fractal dimension of the classical strange rcpcllcr. This 
result has been tested in some systems [7, 8] . As a conse- 
quence, the majority of the eigenfunctions become degen- 
erate with their eigenvalue modulus tending to zero as the 
size of the opening (the number of decay channels) rela- 
tive to % increases. These are the short lived eigenstates, 
which cannot be associated to any classical trapped set 
(instead, they can be related to the trajectories that es- 
cape from the system before the Ehrenfest time). On 
the other hand, the number of remaining eigenstates (the 
long-lived ones) tends to zero. However, they contain the 
most relevant classical information, resembling the clas- 
sical repeller. This was noticed in [9], where they were 
coined quantum fractal eigenstates. Moreover, this in- 
vestigation was recently extended [10] by looking at the 
right and left resonances of the open baker's map. It was 
found that the probability density averaged for several 
right eigenstates is supported by the classical Cantor set 
(the repeller), showing self-similarity both in the q and 
p representation. Finally, in the more specific context of 
optical microcavities, the formation of long-lived scarred 
modes has been observed [3]. This behaviour has been 
associated to avoided resonance crossings. Nevertheless, 
almost nothing else is known about the morphology of 
individual resonances. 

We are interested in the study of quantum systems 



which are classically chaotic. In closed quantum chaotic 
systems, the morphology of the eigenfunctions has been 
extensively studied. One of the most important and 
striking properties is scarring [11]. This consists of the lo- 
calization, i.e., the probability enhancement of given indi- 
vidual eigenfunctions along short unstable periodic orbits 
(POs). This effect has been discovered in the Bunimovich 
stadium billiard [12] and a great amount of work has been 
done since then [13], giving rise to what is known as "scar 
theory" . 

In this letter we explore quantitatively the localization 
properties of resonances. We have studied the overlaps 
of wavefunctions highly localized on the vicinity of POs 
(scar functions) [14, 15, 16] with the eigenfunctions of 
an open quantum system in order to unveil the quantum 
mechanical manifestation of short POs. These values re- 
sult to be higher than when the overlap is calculated with 
the eigenfunctions of the closed system. The U smoothed 
fractal nature does not destroy structures of this kind. 
This effect is even greater when the area of the opening 
grows, thus it cannot be ascribed to a perturbative origin. 
We provide with an interpretation of these results. 

One of the most studied open systems correspond to 
two dimensional torus maps, where a band along the q 
or p directions is cut by means of a projection. The cor- 
responding quantum dynamics is given by a nonunitary 
matrix M = PM (or equivalently M' = MP which is 
related to M by a time reversal operation), where M is 
the closed map and P is the projector on the comple- 
ment of the opening. This quantum evolution is charac- 
terized by decaying eigenstates (pi, whose corresponding 
eigenvalues z\ have complex energies. It is usual to de- 
fine \zi\ 2 = exp (— Ti), where Tj > is called the decay 
rate. We have studied a paradigmatic model of quantum 
chaos, the cat map, which is a linear automorphism on 
the 2-torus generated by the 2x2 symplectic matrix A4, 
modulus 1. We have used 



M = 



2 3 
1 2 



(1) 



When quantizing the torus we have a finite Hilbert space 
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of dimension N — \/2nT% and a discrete N lattice of posi- 
tion and momenta in the unit interval. The quantum cat 
map in the position representation is given by the matrix 
M whose elements are [17] 



M, 



kj 
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exp 
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(2) 



Finally, we choose to apply the projector P after M to 
obtain the nonunitary matrix M that gives the evolution 
of the open cat map. 

The main resource that we use to investigate local- 
ization is the scar function, which not only applies to 
maps but also to general flows. These functions have 
been deeply studied in the literature [14, 15, 16]. They 
are wavcfunctions highly localized on the stable and un- 
stable manifolds of POs, and on the energy given by a 
Bohr-Sommerfcld quantization condition on the trajec- 
tory. We are going to use a formulation suitable for a 
Poincare surface of section, or more directly for maps of 
the 2-torus (examples of this can be found in [15]). We 
define the "tube functions" for maps, \<j\nhc) as a sum °f 
coherent states centered at the fixed points of a given PO 
/i, each one having a phase [16]. Then, a dynamical av- 
erage is performed, and we have the following expression 
for the scar function 
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where T stands for the number of iterations of the map 
up to the Ehrcnfest time Tg = In h/X (A is the Lyapunov 
exponent), and 5 M is the classical action of \x. We have 
used Eq. (3) to construct functions highly localized on 
the vicinity of the periodic points of the closed cat map. 
In Fig. 1(a) we can see the structure of the scar function 
corresponding to the PO given by (q,p) = (0.5, 0.5), one 
of the shortest of this map, for N = 225. The maximum 
probabilities correspond to the darkest regions. Panel (b) 
of the same Figure displays this function in a logarithmic 
scale of gray, showing the way it extends along the stable 
and unstable manifolds of the corresponding orbit. 

In the following we are going to describe the behaviour 
of localization in the open system by means of the max- 
imum overlaps of the scar function with its resonances. 
We explore different values of N and two different shapes 
of the projector P. For simplicity we define a map 
M a = P a M, where P a corresponds to the projection on 
the complement of a band parallel to the p direction of 
width Aq and centered at q = qo . We also define the map 
Mb = PbM, where we now consider two symmetric bands 
each one having a Aq/2 width, and centered at q = qo 
and q = 1 — qo . This is shown in the left and right insets 
of Fig. 3. In Fig. 1(c) and (d) we show the right and 
left eigenstates of the M a map that have maximum over- 
lap with the scar function displayed in Fig. 1(a) (here 



qo = 0.225 and Aq = 0.25). The left resonance localizes 
on the unstable manifold, while the right one does it on 
the stable manifold. The same can be found in panel (e) 
but for the Mb map (go = 0.1625). We can see that the 
symmetric cut localizes the resonance on the stable and 
unstable manifolds of the trajectory. Finally, in panel (f) 
the eigenfunction of the closed cat map with maximum 
overlap with the scar function of panel (a) is shown. 
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FIG. 1: (a) Scar function ofthePO given by (q,p) = (0.5,0.5), 
for N = 225. The horizontal axis corresponds to the position 
q £ [0, 1] and the vertical axis to the momentum p € [0, 1] 
coordinate, (b) Logarithmic version of (a). Black lines corre- 
spond to the stable and unstable manifolds of the orbit, (c) 
Right eigenvector of the M a map which has the maximum 
overlap with the scar function shown in (a), (d) Same as pre- 
vious panel but for the left eigenstate. (e) Same as (c) in the 
Mb map case, (f) Eigenfunction of the closed system having 
maximum overlap with the scar function of panel (a). Sec 
text for more details. 

First, we systematically analyze the behaviour of lo- 
calization as a function of h. For that purpose, in Fig. 
2 we show the maximum overlaps of the scar function 
with the right and left eigenstates of the open cat map, 
as a function of N (for clarity of the exposition we show 
the running average of these values in a window of size 
AiV = 10). It is evident that these values are greater for 
the open system, both for the right (blue dotted) and left 
(red dashed) resonances. On the other hand we can see 
the insets, where the order number ^ max of the eigenstate 
with maximum overlap with the scar function is plotted 
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FIG. 2: (Color online). Maximum overlap x ma x of the scar 
function with the eigenstates of the open cat map as a func- 
tion of N (running average in a window AiV = 10). In the 
insets we show the order number f ma x in ascending eigen- 
value modulus of the right resonance with maximum overlap, 
as a function of N (for the left resonance, results are sim- 
ilar). Panel (a) corresponds to the M a map, while the Mb 
map results are shown in (b). In all panels: the solid black 
lines correspond to the maximum overlap for the closed map, 
and the blue dotted and red dashed lines correspond to the 
maximum overlaps with the right and left resonances of the 
open maps, respectively. 



vs. N (they were ordered in ascending eigenvalue mod- 
ulus). It is clear that the maximum overlap corresponds 
to resonances with the smaller decay rates (larger eigen- 
value moduli). This guarantees that we are looking at 
wavefunctions which have a support on the classical re- 
peller and do not belong to the null subspace. In Fig. 
2(a) we can see the results for the M a map, while in 
(b) they correspond to Mb, where go and Aq values are 
taken the same as in the particular case of Fig. 1. The 
difference between the two maps turns out to be very 
important. In fact, the greater overlaps were obtained 
when opening in a symmetric fashion rather than with 
a single strip. Finally, we mention that the overlaps of 
the scar function with the left or the right resonances of 
the open map differ. These eigenstates are supported by 
different trapped classical sets, so in principle there is no 
reason for them to coincide. Anyway we think that the 
detailed explanation of this difference is an interesting 
open problem. 



But then a natural question arises: how does the rela- 
tionship between the shape and the size of the projection 
influences the intensity of scarring? For instance, this is 
relevant if we want to obtain highly localized resonances 
with the minimum amount of losses. This happens in 
many applications, the cases of two-dimensional billiards 
that can be used as optical microcavities for lasers or that 
can be attached to perfect leads, being some examples. 
To answer this we have further investigated the behaviour 
of localization by fixing the h value, and studying how 
the width of the opening influences it for both, P a and Ph 
operators. The results are shown in Fig. 3, where we dis- 
play the average of x max taken from N — 350 to N = 360, 
as a function of the width of the opening Aq. In all cases 
we take q Q = 0.225 for M a and q Q = 0.1625 for M b . 
The overlaps were calculated with the right (blue dot- 
ted) and left (red dashed) eigenstates. The lower curves 
correspond to M a , while the upper ones to M b . We have 
found that not only the overlap in general increases with 
the size of the opening, but also that this effect is greater 
due to the symmetrization. 
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FIG. 3: (Color online). Maximum overlap x ma x (average from 
N = 350 to N — 360) of the scar function with the right 
(blue dotted) and left (red dashed) eigenstates of the open 
cat map as a function of the size of the opening Aq. The 
lower curves correspond to M a , while the upper ones to M&. 
The solid horizontal line stands for the value corresponding 
to the closed cat map. Left and right insets illustrate the 
projectors P a and Pt- 

But this seemingly greater scarring effect in open sys- 
tems should be interpreted in the proper context. In 
order to do this we will analyze the weight that these 
long-lived resonances have in the whole spectrum, and 
relate it with typical time scales of the system. This is 
given by a connection between the fractal Weyl law and 
the Ehrcnfest time T a = In (O) / A (with O the number of 
open channels, and A = 1.31 in our case), first obtained 
in [7] . There, it was found that the fraction of resonances 
with decay rate T smaller than a fixed value r* < 1/T Q 
behaves like N 7 /N = exp(-T /T d ) = 0(-V(^)) ) where 
T d = N/O is the so-called "dwell time" (T d large). We 
have numerically confirmed the validity of this predic- 
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tion in our system by fitting the data with the expression 
Nry/N — aN~ b . We show three cases in Fig. 4, where 
we have taken Fy = 0.71 in all of them. The upper curve 
corresponds to M a with an opening defined by q Q = 0.125 
and Aq — 0.05, showing a fitted bf = 0.032 that agrees 
with the theoretical b t h — 0.038. The middle curve corre- 
sponds to go = 0.225 and Aq = 0.25 with bf = 0.181, and 
the lower one corresponds to with go = 0.1625 and the 
same Aq with bf = 0.2, being bth — 0.191 for both. In all 
cases this fraction goes to zero, leaving a small amount 
of classically meaningful eigenstates. This directly im- 
plies a persistent localization effect on the few remaining 
"fractal eigenfunctions" . However, these greater overlaps 
could correspond to a normalization difference with re- 
spect to the closed system. In fact, the effective size of 
the Hilbert space of the open system is smaller. Then, 
if we can go further and claim that these results mean a 
true enhancement of scarring, remains to be determined 
[19]. 
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FIG. 4: (Color online). Fraction of eigenstates N-y/N whose 
decay rate Y is smaller than Tj = 0.71, as a function of 
N. The red lines corresponds to the theoretical expression 
Q(~ 1 /(^ T d)) (•ggg text for details). The upper curve corre- 
sponds to M a with an opening defined by go = 0.125 and 
Ag = 0.05, the middle one by go = 0.225 and Ag = 0.25, and 
finally the lower one corresponds to Mb with go = 0.1625 and 
Ag = 0.25. 

In summary, we have found that there is a greater over- 
lap of the scar functions with the resonances of an open 
system compared to the closed one. The fractal struc- 
ture of the eigenstates has been widely studied, motivat- 
ing their denomination as "quantum fractal eigenstates" . 
However, this significant alteration of the morphology of 
the eigenfunctions with respect to the analog closed sys- 
tem does not destroy the localization around POs. We 
think that this is due to the fact that the pruning of orbits 
that escape through the openings before the Ehrenfest 
time leaves parts of the stable and unstable manifolds. 



These remaining parts are enough to support the quan- 
tum structure associated to the scar function. However, 
the way they interfere in order to construct the same 
object than the smooth manifolds of the closed system 
remains unknown. Also, the scarring grows with the size 
of the opening ruling out any perturbative explanation 
for this. Moreover, this phenomenon is persistent, in the 
sense that it survives in the vanishing fraction of long- 
lived resonances as N grows. In future studies [19] we 
will focus on these open questions. 
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